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INTRODUCTION 


An  artillery  round  with  a  combustible  cartridge  case  is  fired 
from  a  weapon  using  a  control  system  which  loads  the  round,  aims  the 
weapon  and  fires  the  weapon,  all  automatically.  At  least  two  rounds 
of  ammunition  are  fired  in  this  fashion  and  some  smoldering  residue 
from  the  preceding  round  may  remain  in  the  chamber  of  the  weapon  when 
a  round  is  loaded.  Let  R  be  the  conditional  probability  that  the 
sidewall  of  the  cartridge  case  of  the  chambered  round  is  not  burned 
through  prior  to  firing,  given  that  it  is  instantaneously  ignited  by 
smoldering  residue.  The  problem  is  to  find  both  a  point  estimate  and 
a  95%  lower  confidence  limit  for  R,  using  information  concerning  the 
gun  cycle  time  and  data  on  cartridge  case  burn-through  time  obtained 
from  laboratory  tests. 

There  was  sufficient  gun  cycle  time  data  available  to  justify  the 
assumption  that  the  elapsed  time  between  chambering  and  firing  a  round 
is  normally  distributed  with  (true)  mean  y2  =  2.9  seconds  and  standard 
deviation  a2  =  0.13  seconds. 

An  experiment  was  conducted  to  estimate  the  statistical  distribution 
of  cartridge  case  sidewall  burn-through  times.  One  hundred  and  fifty 
samples  cf  sidewall  material  were  taken  from  several  cartridge  cases 
and  tested.  Each  sample  of  sidewall  material  was  ignited  and  the 
elapsed  time  between  ignition  and  burn-through  was  measured  by  three 
observers  using  stop  watches.  No  data  were  obtained  for  two  samples 
and  there  were  some  missing  data  for  some  of  the  remaining  samples. 
Unfortunately,  because  of  the  developmental  nature  of  the  item  under 
consideration,  it  was  necessary  to  sample  cartridge  cases  from  only 
one  lot  and  assume  that  this  sample  was  randomly  selected  from  the 
conceptual  population  consisting  of  all  such  cartridge  cases  which  will 
be  manufactured.  It  was  suggested  that  the  validity  of  this  assumption 
should  be  verified  by  further  testing  when  a  sample  which  is  more 
representative  of  the  production  item  can  be  selected. 


PRELIMINARY  STATISTICAL  ANALYSIS 

The  burn-throuqh  data  were  analyzed  assuming  a  one  way 
classification,  components  of  variance  (random  effects),  analysis  of 
variance  model  with  unequal  number  of  observations  per  cell.  The 
component  of  variance  attributable  to  observers  was  much  smaller  than 
the  among  samples  component  of  variance  (0.08  sec.2  vs.  2.23  sec.2) 
and  was  not  statistically  significant  at  the  a  =  0.001  level  o-f  significance. 
It  was  concluded  that  the  precision  of  measurement  resulting  from  the 
use  of  observers  with  stop  watches  was  adequate. 

The  data  were  used  to  test  the  hypothesis  of  normality  of  the 
distribution  of  burn-through  times,  a  requirement  of  the  Church-Harris 
procedure  used  for  estimating  R.  A  chi-square  qoodness  of  fit  test 
rejected  this  hypothesis  at  the  0.05  level  of  significance  but 
accepted  it  at  the  0.01  level.  Since  the  chi-square  test,  which 
is  relatively  insensitive  to  departures  from  normality  in  the  region 
of  the  tails  of  a  distribution  was  inconclusive,  the  statistic 
bj  =  n[z(x1-  -  x)3]2[z(xi  -  x)2]"  was  calculated  and  used  to  test 
the  hypothesis  that  the  distribution  is  not  skewed.  This  hypothesis 
was  rejected  at  the  0.02  level  of  significance  (a  two  tail  test  with 
0.01  probability  in  each  tail  was  used)  so  it  was  inferred  that  the 
distribution  is  positively  skewed.  Next,  Craig's  procedure  [2]  was 
used  to  determine  which  member  of  the  Pearson  system  of  frequency  curves 
best  describes  the  data.  It  was  found  that  the  Pearson  Type  III  curve 
(a  gamma  density  function)  fits  the  data  best.  From  Carver's  table  of 
the  standardized  Type  III  function  [3],  it  was  verified  that  the  lower 
tail  of  the  Pearson  Type  III  curve  contains  less  area  in  the  interval 
-°°  <  x  <_  v  +  a  than  a  normal  curve  with  the  same  mean  and  variance. 

This  indicated  that  a  normality  assumption  would  leave  to  conservative 
point  and  interval  estimates  of  R,  i.  e.,  if  the  estimates  are  biased, 
the  bias  will  be  such  that  R  is  underestimated. 
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For  the  purpose  of  estimating  R,  It  was  assumed  that  cartridge 
case  sidewall  burn- through  time  Is  normally  distributed  with  estimates 
of  the  mean  and  standard  deviation  of  the  distribution  being  X  =  9.71 
seconds  and.  S  =  1.49  seconds,  respectively. 

APPLICATION  OF  A  MODIFIED  CHURCH-HARRIS  PROCEDURE 

The  procedure  used  for  estimating  R  Is  based  on  the  work  of  Church 
and  Harris  [1].  Let  the  random  variable  X  be  the  cartridge  case  sidewall 
burn-through  time  In  seconds  and  the  random  variable  Y  be  the  gun  cycle 
time  in  seconds.  Assume  that  X  and  Y  are  statistically  independent  and 
both  normally  distributed.  Introducing  the  notation 


E(X)  =  V1 
VAR(X)  =  aj 
E(Y)  =  y2 


VAR(Y)  = 

and  defining  the  random  variable  W  =  Y  -  X,  it  follows  that  W  is  normally 
u.stributed  with  E(W)  =  U2  -  Pl  and  VAR(W)  =  +  a*.  Then 


R  =  P{Y  <  X)  =  P{Y  -  X  <  0}  =  P{W  <  0}. 
We  next  make  the  transformation  Z  = 


W  -  (u2  -  ^i) 


so  that  Z  is 


distributed  normally  with  E{Z)  =  0  and  VAR(Z) 
R  =  P{W  <  0}  =  P  |  Z/^2  +  ^2  +  Vz  ~  Vi  <  oj 


(1) 
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where  *(•)  Is  the  standard  normal  cumulative  distribution  function. 
Substituting  the  known  values  w2  and  c2  and  the  estimated  values 
uj  =  X  and  a2  a  s2  Into  (1)  yields  the  point  estimate 


R  = 


X  -  la  2 
W  +  a\ 


) 


Now,  having  a  point  estimate  of  R,  we  proceed  to  approximate 
the  probability  distribution  of  the  random  variable  R  and  use  this  to 
construct  an  approximate  100(1  -  y)2S  lower  confidence  limit  for  R. 

In  doing  so,  we  make  use  of  the  fact  that  9  is  normally  distributed 
with  E(X)  =  vi  and  VAR(X)  =  of/n;  S2  is  asymptotically  normally 
distributed  with  E(S2)  =  a2  and  VAR(S2)=2c^/(n-l ) ;  and  that  9  and  S;' 

are  statistically  independent. 

20i 

Let  I  =  S2  -  of  so  that  E(T)  =  0  and  VAR(T)  =  — j- 


and  define 

X  —  p 2  X  -  P2 

V  =  -  = . . - 

/S2  +  a2  /a 2  +  o|  +  T 


Expanding  V  in  a  Taylor's  series  about  the  point  [E(X),  E(T)]  =  (ui »  0) 
we  obtain 

9  -  v2  1  (v\  ~  V2)T  t 

v  =  - -  i  t  0(-h  . 

/a i  +  °2  ^  +  °2) 

with  probability  one. 

Because  of  the  indenendence  of  9  and  T,  the  distribution  of  V  is 
asyasymptotically  normal  with 

E  ( V)  =  - - 
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and 


VAR(V)  * 


°?(ui  -  uz)  "I 
(n-l)(of  +  af)2j 


Using  S2  to  estimate  a2  and  X  to  estimate  v\  we  obtain  as  an  estimate 


of  the  standard  deviation  of  V 


a 


V 


$ _ 

^  +  o22 


S2(K  -  y2)2  I 
(n-l)(S2  +  a2 )2  J 


and 


Since  R  = 


=  *[E(V)] 


it  follows  that 

P  |r  >  *[V  -  4>  (1  -  y)ov]|  «  1  -  y 

We  use  (2)  to  obtain  the  1 00 (  1  -  y)  %  lower  confidence  limit  for  R  as 


LCL1_y=$[V  -  «  (1  -  Y)ay]. 
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COMPUTATIONAL  RESULTS 


For  the  combustible  cartridge  case  problem  we  want  a  95%  lower 
confidence  limit  for  R.  Since  <i>“ 1  ( . 95 )  =  1.645,  this  limit  is  qiven 
by 

LCL  95  =  *(V  -  1.645ov). 

A  computer  program  (Appendix  A)  was  prepared  to  calculate  V,  R,  ay 

and  a  95%  lower  confidence  limit  for  R.  In  this  particular  problem  a 
95%  lower  confidence  limit  of  1  -  10"4  was  considered  satisfactory  to 
assure  the  margin  of  safety  reauired.  The  95%  lower  confidence  limit 
was  determined  to  be  1  -  0.2  X  1C*4. 

Next,  the  question  of  how  much  the  mean  gun  cycle  time  can  be 
increased  and  still  leave  a  95%  lower  confidence  limit  of  1  -  1CT4 

for  R  was  considered.  To  answer  this  question  it  was  assumed  that  the 

coefficient  of  variation  of  the  random  variable  Y,  02^2*  remained 
constant  as  v2  increased,  from  2.9  seconds  to  3.9  seconds,  in  steps  of 
0.1  seconds.  V,  R,  and  the  95%  lower  confidence  limit  for  R  were 
calculated  at  each  step.  It  was  fv,u..d  that  the  mean  gun  cycle  time 

can  be  increased  as  much  as  0.5  seconds  without  the  95%  lower  confidence 

limit  for  R  falling  below  1  -  10"4. 

SUMMARY  AND  CONCLUSIONS 

The  application  of  the  Church-Harris  technique  to  the  combustible 
cartridge  case  Droblem  can  be  summarized  as  follows:  The  analysis 
depended  on  two  critical  assumptions;  a)  the  sample  of  cartridge  case 
sidewall  material  used  to  obtain  burn-through  time  measurements  was  a 
random  sample  from  the  conceptual  population  consisting  of  all  cartridqe 
cases  of  the  same  type  which  will  be  manufactured  in  the  future  and  b) 
the  distributions  of  X  and  Y  are  normal.  Assumption  a)  seemed 
questionable  and  it  was  suggested  that  further  testing  be  done  to 
verify  it.  Assumption  b)  was  not  satisfactorily  established  by  the 
data  but  the  analysis  indicated  that,  if  the  assumption  is  not  valid 
the  inferences  drawn  from  the  study  will  be  on  the  safe  side,  i.  e., 
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R  win  be  underestimated.  The  conditional  probability  that  the 
sidewall  of  a  cartridge  case  of  a  chambered  round  will  not  be  burned 
through  prior  to  firing,  given  that  It  Is  Instantaneously  Ignited  by 
smoldering  residue  remaining  from  a  round  previously  fired,  was  estimated 
to  be  0.9999971  with  a  95*  lower  confidence  limit  of  0.9999779.  It 
was  also  determined  that  the  mean  gun  cycle  time  can  be  Increased  as 
much  as  0.5  seconds  (a  17*  Increase)  without  the  95*  lower  confidence 
limit  exceeding  0.9999,  provided  that  the  coefficient  of  variation 
remains  constant. 


REFERENCES 


[1]  Church,  J.  D.  and  Harris,  B.  (1970).  "The  Estimation  of 
Reliability  from  Stress-Strength  Relationships".  Teahnometrias  12, 
49-54. 

[2]  Craig,  C.  C.  (1936).  "A  New  Exposition  and  Chart  for  the  Pearson 
System  of  Frequency  Curves".  The  Armais  of  Mathematical  Statistics ,  7, 
16-28. 

[3]  Carver,  H.  C.  (1940).  "Statistical  Tables".  Eduards  Brothers,  Inc., 
Arm  Arbor,  Michigan. 


14 


APPENDIX  A 


The  following  is  a  listing  of  a  subroutine  CHAR  (CHurch-HARris) 
and  a  representative  driving  program  written  for  this  study.  Notice 
that  CHAR  requires  as  Input  X(a  vector  of  data),  y2»  °2»  y  »  N  and 
outputs  X,  S,  R,  ov,  and  LCL.  Subroutine  CHAR  calls  subroutines  FND  and 
FINVND  to  evaluate  the  normal  distribution  function  and  the  inverse  of 
the  normal  distribution  function,  respectively. 

The  output  of  the  sample  program  is  formatted  as  appears  below. 

The  dimension  statement  (DIMENSION  X(500))  appearino  in  the  drivinq  program 
may  be  modified  to  accommodate  the  data  available. 
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o  o  o  non  o  o  o  n  c » rv 


REAL  LCL, MU2 
DIMENSION  X ( 500) 

READ  INPUT  PARAMETERS 

READ! 5« 101)MU2,SIGMA2,GAMMAfN 

101  FORMAT! 3F10. 3, 110) 

READ  OATA 

READ! 5, 102)1X11), I*1,N) 

102  FORMAT! 10X, 3F10. 5) 

CALL  CHAR ! X, MU2,*SIGMA2, GAMMA, N.XBAR, S, RHAT, V,S IGVHAT ,LCL ) 

RHAT  * 1.0-RHAT 
LCL-l.O-LCL 

WRITE  COLUMN  HEADINGS 

WRITE16.103) 

103  FORMAT! 6X, 3HMU2, 7X.6HSI GMA2.6X, 1HN,6X,4HXB AR , 10X , 1HS ,8X,6H1-RHAT , 
19X , 1HV, 8X, 7HS I GVHAT , 6X,  5HGAMMA , 7X , 5H1-LCL > 

MR ITE 1 6* 104) 

104  FORMAT! 1H0) 

WRITE  OUTPUT 

WRITE! 6, 105 )MU2,SIGMA2, N.XBAR ,S,RHAT .V.SIGVHAT .GAMMA, LCL 

105  FORMAT! 2E12.3, I6.7E12.3) 

STOP 


END 


SUBROUTINE  CHAR ! X, MU2, S IGMA2 .GAMMA ,N, XBAR , S , RHAT , V , S IGVHAT ,  LCL ) 

CHAR 

1 

REAL  LCL.MU2 

CHAR 

2 

DIMENSION  XII) 

CHAR 

3 

XN»N 

CHAR 

4 

XBAR-0 

CHAR 

5 

SX2«0 

CHAR 

6 

DO  11  1*1, N 

CHAR 

7 

X  BAR'XB AR  +X ( I ) 

CHAR 

8 

SX2*SX2+X!I)*X(I) 

CHAR 

9 

11  CONTINUE 

CHAR 

10 

XBAR*XBAR/XN 

CHAR 

11 

S2*! SX2-XN*XBAR**2 )/!XN-l.O) 

CHAR 

12 

S*SQRT!S2) 

CHAR 

13 

VAR=S2*SIGMA24*2 

CHAR 

14 

V* ( XBAR-MU2 ) /SORT! VAR ) 

CHAR 

15 

RH*T=FND!V) 

CHAR 

16 

S IGVHAT  =  1 S/SQRT ( VAR ))*SQRT(1.0/XN*0.5*S2*(XBAR-MU2)**2/( (XN-1.0) 

CHAR 

17 

14VAR**2 ) ) 

CHAR 

18 

LCL*FND(V-FINVN01 1,0-GAMMA )*SIGVHAT) 

CHAR 

19 

RETURN 

CHAR 

20 

END 

CHAR 

21 

16 

_ 


